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We have measured, in real time, individual quasiparticle tunneling in a single Cooper pair tran- 
sistor, using rf reflectometry on the supercurrent branch. We have extracted the even-to-odd and 
odd-to-even transition rates directly by analyzing the distributions of dwell times in the even and 
odd states of the transistor. We discuss both the measurement and analysis techniques and report 
on the temperature and gate dependence of the quasiparticle tunneling rates. 
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Superconducting single-charge devices form an impor- 
tant subclass of quantum circuits, of which perhaps the 
best known is the Cooper-pair box (CPB) qubit.l^ In 
the CPB and several similar circuits, uncontrolled quasi- 
particle (QP) tunneling in the system can degrade their 
performance, causing the operating points of these de- 
vices to shift stochastically on time scales comparable 
to the measurement times. In CPB qubits, this can 
have the ultimate effect of reduced visibility of coher- 
ent oscillations, and may become an important decoher- 
ence mechanism as coherence times improve In a re- 
lated circuit, the single Cooper-pair transistor (SCPT), 
this uncontrolled 'poisoning' of the supercurrent can be 
a limiting factor in its use as an electrometer, while QP 
tunneling in Cooper-pair pumps^ severely inhibits their 
applicability as metrological current sources. While this 
problem has already been studied in some detail with 
other methods, in this paper we present rf measurements 
of SCPTs that allow us unprecedented access to the dy- 
namics of QP tunneling in this class of circuits. As shown 
below, we have been able to monitor this process in the 
time domain and, by analyzing the resulting telegraph 
signal, determine the quasiparticle tunneling rates in a 
straightforward manner. 

Quasiparticle poisoning has been extensively investi- 
gated in the past both theoretically and experimentally 
in the SCPT.i2i2ii2iiiii2iiii4 In those experiments the 
relative population of the odd (an excess QP on the is- 
land) and even (no excess QP on the island) states of the 
transistor was inferred, either by measuring the average 
charge on the island or by measuring the distribution 
of the current /gw at which the transistor switches out 
from its supercurrent branch. These experiments how- 
ever do not measure the even-odd and odd-even transi- 
tion rates, 7eo and 7oe directly, whereas these rates are of 
particular interest in the context of charge qubit decoher- 
ence. We note that switching current measurements can 
in principle provide direct information about QP tunnel- 
ing rates^ but are sensitive only to QP tunneling events 
that reduce the critical current of the transistor: at gate 
charges Ug = CgVg/e £ [0.5, 1.5] these are even^ odd 
transitions {Cg and Vg are, respectively, the gate ca- 



odd— !■ even transitions. Furthermore, inference of tun- 
neling rates from these measurements is complicated by 
the evolution of the system during the current ramp. 

Here we present the first direct measurement of QP 
tunneling rates in a zero-biased SCPT using time-domain 
rf reflectometry with microsecond temporal resolution. 
In contrast to previous experiments, our measurement is 
sensitive to both even— > odd and odd-^- even transitions. 
By measuring the impedance of the transistor biased on 
its supercurrent branch, we can read out the charge on its 
island without switching to the voltage state, thus avoid- 
ing additional pair breaking in the device."'^^ Reading out 
the charge in this manner also avoids the added compli- 
cation of re-initializing the system after each switching 
event. 

To perform these experiments we embedded the SCPT 
in a resonant circuit, whose resonance frequency is de- 
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FIG. 1: (Color online) A schematic of the experimental setup. 
RF power is fed to the circuit via a directional coupler after 
passing several stages of attenuation and low-pass (LP) filter- 
ing. The refiected power is amplified at 4 K, with additional 
amplifiers at room temperature. Cc = 3.3 pF, La = 18 nH, 
and Cs = 6.2 pF. A 22 pF capacitor in series with Ls, and an 
rf choke Lc = 820 nH allow dc biasing of the transistor. Inset 
(a): |rp vs. frequency for device H, at difi'erent gate charges. 
Inset (b): |rp vs. rig at 515 MHz. The power incident on the 
resonator was Pin 



10"'^ W (-120 dBm). 
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termined by the parallel combination of Cg, Lg (Fig.^, 
and the source-drain impedance of the transistor. The 
impedance of the zero-biased SCPT includes the Joseph- 
son inductance)i& Lj(ng) ~ h/2elc(ng), where Icing) is 
the charge dependent critical current. We found that 
this impedance additionally contains a real part, Rj{ng), 
likely due to fluctuations in the Josephson phase induced 
by the device's finite temperature and the local environ- 
mental impedanceiiiiSiiS An rf carrier signal probes the 
circuit near its charge-dependent resonance frequency so 
that the power reflected from the resonator is a strong 
function of the island charge. A measurement of the re- 
flection coefficient |rp — Pout/^in is therefore a mea- 
surement of the charge on the island, Fig. ^a), and its 
temporal evolution measures the charge dynamics in the 
transistor. The circuit parameters were chosen to set 
its resonance near 500 MHz, the center of our circulator 
and amplifier band, and the probe power Pin was set low 
enough to keep Lj in the linear regime. Also note that 
the design of the circuit allows us also to characterize the 
device by use of conventional dc techniques. 

Here we show data from two devices fabricated on 
the same chip^SS using double-angle Al deposition. By 
oxygen-doping the Al in the first deposition step, we 
have produced films with different superconducting gaps, 
Al = 225 ^eV (20 nm thick) and A2 = 190 ^eV (40 nm 
thick). In the first SCPT, which we will refer to as de- 
vice H, the film with the higher gap Ai formed the island 
of the transistor, and the leads had the smaller gap A2. 
In the second SCPT, which we call device L, the gap of 
the island (A2) was lower than that in the leads (Ai).-^ 
Otherwise the two devices were nearly identical: for de- 
vice H, Ec = 170 fieV, i?N = 21.2 kl7, and for device L, 
Ec — 177/ieV, i?N = 21.3MI. Despite the similarity of the 
two devices, their different gap profiles lead to strikingly 
different QP 'poisoning' dynamics — device L tends to 
trap quasiparticles on its island far more effectively than 
device Hi^ 

The eigenenergies of the transistor's Hamiltonian, and 
correspondingly its critical current, are 2e periodic in 
the external charge applied to the island. Quasiparticle 
tunneling onto and off of the island shifts the net island 
charge by one electron, offsetting the modulation charac- 
teristic by half a period. Measurements of the transistor 
on time scales greater than the characteristic poisoning 
time will consequently show an average of its properties 
at Ug and rig ± 1, weighted respectively by the even and 
odd state probabilities. Faster measurements will show 
transitions between these charge parity states. 

In device H we observed no QP tunneling events and 
it appears 'poisoning free'. However, from the device pa- 
rameters and the QP trapping model of Ref. 0, we do 
expect fast QP transitions to occur near Ug = 1 and con- 
clude that their rates fall outside the bandwidth of our 
measurement. Our maximum measurement bandwidth^^ 
of 3 MHz sets an upper bound on the lifetime of the odd 
parity state of the transistor, Tq < 0.3 fj,s. In addition, 
we observe no significant averaging of the reflected power 
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FIG. 2: (Color online) (a) Intensity plot of the reflection 
coefiicient |r|^ (dB) in device H, measured as a function of 
gate charge and frequency at T = 25 mK (Pin — 10"^^ W). 
(b) Calculated |rp with Lj{ng) and Rj{ng) obtained from 
fits, (c) Fitted values of Lj as a function of gate (circles), 
and the calculated Lj(ng) (line), (d) Fitted values of Rj. 



away from its expected 2e modulation [Fig. QJb)], sug- 
gesting that the transistor is predominantly in the even 
parity state. This requires that the even state lifetime 
be greater than Tq by at least an order of magnitude. 

To > 3 /iS. 

We used the apparently 'clean' 2e-periodicity of device 
H to verify the operation of our rf circuit. Figure [21[a) 
shows the measured reflection coefficient jFp as a func- 
tion of frequency and gate charge in device H. At each 
gate charge, we fitted the experimental |r(tt')p with our 
circuit model, allowing Lj and Rj to vary while holding 
the other circuit elements fixed close to their nominal 
values. The fitted response is shown in Fig. |2Ib). Fig- 
ure [2Ic) shows Ljijig) obtained from the fits (circles), 
overlaid with the expected charge dependence of Joseph- 
son inductance in this device, calculated as in Ref. hj 
(line). We found that the charge-dependent loss, Rj{ng), 
modulates between approximately 20-200 fl [Fig. |2Id)], 
values that are somewhat lower than the phase diffu- 
sion resistance we measured at dcm^ For the purpose 
of the present experiments, our simple model for the 
SCPT's impedance — an inductor in series with a re- 
sistor — agrees with the data sufficiently well; a full de- 
scription of this impedance at arbitrary temperature, fre- 
quency, ac drive, and dc bias will require further exper- 
imental studyiiS As an electrometer, the device's best 
charge sensitivity was 5q = 5.2 x 10~^ c/a/Hz, as es- 
timated from the sideband signal-to-noise ratio with a 
small amplitude signal (0.025 e rms) at / = 1.5 MHz 
applied to its gateJ^ These numbers are similar to the 
performance of conventional (dc-biased) rf-SET's,^^ and 
are likely to improve with careful design of the rf circuit 
and the use of lower noise cryogenic amplifiers. 

We now turn to device L, where the gap profile was 
reversed. Since the superconducting gap of the island 
in this device was smaller than that of the leads, the 
island trapped quasiparticles for longer periods of time 
than in device H, making the QP capture and emission 



3 




35 
30 
25 
_ 20 
15 
10 
5 




n=1.0 


n 


0.85 




• 




■ ^ 


V 





25 



75 



125 
T(mK) 



175 



FIG. 3: (Color online) (a) |rj vs. time measured in device 
L at rig = 1. The resonator was probed at / = 505 MHz 
with Pin = 1.6 X lO"^'^ W (-118 dBm). (b) A histogram of 
the time record. The peaks correspond to the even and odd 
states of the SCPT. (c) The distribution of observed lifetimes 
in the even and odd states at = 1. (d) Intensity plot of 
time record histograms as in (b), as a function of Ug. Darker 
shades represents higher count number. Note the threshold 
for QP tunneling at n^'" 0.35 (dashed). 



process slow enough to be observed in the time domain. 
Figure |3{a) shows |rp as a function of time, recorded 
at rig — 1. The reflected power switches randomly be- 
tween two levels, seen clearly as distinct peaks in a his- 
togram of the time record. Fig. ^h), where the higher 
level corresponds to the odd state of the transistor with 
a QP trapped on the island, and the lower one corre- 
sponds to the even state. The typical time record in this 
experiment was one second long (2 x 10^ samples), en- 
compassing well over 10"* tunneling events in each trace. 
In Fig. IS^d) we show the gate dependence of the time 
record histograms, revealing not only the rf response of 
the resonator [compare to Fig.^b)], but also the relative 
populations of the even and odd states as a function of 
Ug. Quasiparticle transitions were observed in a range 
of gate charges for which the even-odd energy difference 
SEeo = SEch + SA > 0, where SEch is the difference 
in charging energy and SA = 35 fieV is the gap differ- 
ence between the leads and the island. The calculated 
threshold n^'' ~ 0.35, for which SE^o = 0, is shown in 
Fig.Hd). 

We determined the QP tunneling rates from the dis- 
tribution of dwell times in each of the states. Fig. El^c). 
These distributions are found by operating on the mea- 
sured telegraph signal [Fig.|3^a)] with a change-point de- 
tection algorithm2Si2L2& and histogramming the time in- 
tervals between detected QP transitions. However, the 
experimentally observed transition rates, 7*^ and 7*^, 
which are obtained from exponential fits to the data in 
Fig. |3^c), underestimate the true rates in the underlying 
system, and must be corrected for the finite bandwidth 
of our measurement QP tunneling events whose du- 
ration is shorter than the characteristic resolution time 



FIG. 4: (Color online) Tc (squares, triangles) and To (dia- 
monds, circles) as a function of temperature, at = 1 and 
Ug — 0.85. The data shown was corrected for finite band- 
width effects using the inverted Eq. (0. The solid (dashed) 
line is a fit to the Ug = 1 {rig = 0.85) data, incorporating the 
effects of electron heating by the probe rf signal, giving an 
activation barrier of 112 /xeV (48 /jeV). 



of our receiver cannot be detected, causing the lifetime 
histograms to be artificially skewed towards longer times 
(see also Ref. 's^. We have shown previously that in 
time-domain measurements with a finite-bandwidth re- 
ceiver the observed rates are given by22, 



7o. 



^ (a - -47r7eo(oc)) 



(1) 



where A = 700 + 7oc + 7r , the bandwidth of the receiver 
is 7r = T"-'^, and 7eo and 700 are the true even-odd and 
odd-even rates in the underlying system. For the data 
shown here, = 5 ± 1 ^s. 

In Fig. ^ we show the temperature dependence of the 
even and odd state lifetimes of the transistor, ~ 
and To = ^ from which we can deduce the energetics 
of the poisoning process. The values shown in this figure 
were obtained from the observed rates by inverting Eq. 

to correct for finite bandwidth effects. The error bars 
are dominated by the uncertainty (20 %) in t^. To verify 
the transformation in Eq. we varied the bandwidth of 
our measurement and extrapolated the observed lifetimes 
to T]. — > (infinite receiver bandwidth). Doing so, we 
obtained Tc = 3.0 ± 0.5 ^s and Tq = 26.0 ± 0.7 fxs at 
rig = 1 and T = 25 mK, in agreement with the data in 
Fig.H 

From Fig. 01 we see that Tq, the even-state lifetime, 
does not change significantly with temperature. This 
temperature independence implies that the even— > odd 
transition does not require activation, confirming the ex- 
istence of a population of nonequilibrium quasiparticles 
on the leadsii^ In contrast, Tq changes dramatically with 
temperature above ~ 80 mK, and is clearly activated. 
At lower temperatures, where the electron-phonon ther- 
mal coupling is poor, the power dissipated in the device 
during the measurement causes the electron temperature 
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and the odd-state lifetime to saturate. We fitted the odd 
state hfetime to Tq oc exp{AE/kBTi,), where is the 
electron temperature given by = (P/T,V + T^) .^^ 
Here P is the power dissipated in the transistor, S = 
2 X 10^ W/m3K^ F - 3.2 X lO^^i ^^3 jg jgi^j^^j ygj, 

umc, and T is the temperature at the mixing chamber 
of our dilution refrigerator. The best fit at = 1 (solid 
line in Fig. Q is obtained with P = 1.6 x 10~^^ W, cor- 
responding roughly to the incident rf power, and with an 
activation barrier AE = 112 ^eV. For Ug = 0.85, we find 
AE = 48 /ieV, and the fit is obtained with a lower power 
P = 2.4 X lO"""^^ W, consistent with a higher reflection 
coefficient of the circuit at this gate. 

If the odd^ even transition is due to the odd QP leav- 
ing the island, we would expect the activation barriers 
to be higher than our experimental values. We esti- 
mate that in our devices AE = SE^o — 160 /xeV and 
140 /ieV at Ug — 1 and 0.85, respectively."^^ These esti- 
mates, however, do not take into account quantum cor- 
rections to the even-odd energy difference due to virtual 
QP transitionsSS*^ or higher-order (cotunneling) pro- 
cesses. For example, a second order process involving 
the simultaneous tunneling of a second QP into, and a 
Cooper pair out of the island may dominate the odd-even 
rates. We expect these corrections to be important in our 
devices, whose tunnel resistance i?N < h/e'^ puts them 
in the strong tunneling regime. We have also observed 
activation barriers that are smaller by nearly a factor of 
two than those expected from the orthodox theory, in 
switching current measurements in other SCPT devices 



having i?N = 18 kfl and Ec — 115 /leV. 

To conclude, we have measured single-Cooper-pair 
transistors using time-domain rf reflectometry on their 
Josephson branch. In one of the devices, whose island's 
superconducting gap was higher than that of its leads, 
we observed clean 2e periodicity of the reflected rf power 
as a function of gate charge. Quasiparticle transitions in 
this device, if they occurred, were too fast to be detected. 
In the other device, whose island had a smaller gap but 
was otherwise nearly identical to the first, we measured 
individual QP tunneling events in real time. We found 
the distribution of dwell times in the even and odd states 
of the transistor, and from it extracted the QP tunnel- 
ing rates. This represents the first direct measurement 
of single-electron tunneling rates in a metallic transistor; 
these rates set a lower bound on the decoherence rates of 
superconducting charge qubits.^ The temperature inde- 
pendence of 7oo near rig — 1 indicates that the odd state 
of the transistor is energetically favorable. This obser- 
vation, and the value of rig'', the gate charge at which 
the odd state population becomes significant, are consis- 
tent with the picture of non-equilibrium QP population 
residing on the leads (Ref. 0). The temperature depen- 
dence of the odd—!- even transition rates, however, yield 
activation barriers that are smaller than expected, possi- 
bly due to quantum corrections to the even-odd energy 
difference or to higher-order tunneling processes. 

We thank M.A. Sillanpaa and K.W. Lehnert for helpful 
comments. 
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